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Accuracy of throwing in games and sports is governed by how errors at projectile 
release are propagated by flight dynamics. To address the question of what governs 
the choice of throwing strategy, we use a simple model of throwing with an arm 
modelled as a hinged bar of fixed length that can release a projectile at any angle and 
angular velocity. We show that the amplification of deviations in launch parameters 
from a one parameter family of solution curves is quantified by the largest singular 
value of an appropriate Jacobian. This allows us to predict a preferred throwing 
style in terms of this singular value, which itself depends on target location and the 
target shape. Our analysis also allows us to characterize the trade-off between speed 
and accuracy despite not including any effects of signal-dependent noise. Using 
nonlinear calculations for propagating finite input-noise, we find that an underarm 
throw to a target leads to an undershoot, but an overarm throw does not. Finally, 
we consider the limit of the arm-length vanishing, i.e. shooting a projectile, and 
find that the most accurate shooting angle bifurcates as the ratio of the relative 
noisiness of the initial conditions crosses a threshold. 
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1. Introduction 

What governs the selection of a style when throwing a ball into a bin, say overarm 
versus underarm? At first the question seems ill-posed since of the arbitrarily large 
number of possible choices and variables: people prefer overarm over underarm or 
vice versa depending on age, gender, culture, training etc. However, the question 
becomes better posed if we consider just the physics of projectile release with a 
prescribed velocity and launch angle to reach the target. This is because there is a 
one parameter family of solutions; for almost any given launch angle, we can pre- 
scribe a launch velocity. How then might we choose from this family of trajectories 
? Since any errors introduced in the initial conditions are not uniformly amplified 
by all trajectories, a natural possibility is that we ought to choose that trajectory 
that least amplifies any initial errors. This leads us naturally to the subject of error 
propagation and amplification in dynamical systems with antecedents that go back 



a century or more flPoincare 1912 Hopf 1934) but continue to have implications 



for prediction and decision making. Here, we consider the simple task of throwing 
for two reasons. First, the task decouples the internal (neural) decision making from 
the physics of projectiles, i.e. the task decouples planning from feedback control. In 
the absence of wind, there is no noise during the projectile flight, only amplification 

f mvenkadesan@gmail.com 
\ lm@seas.harvard.edu 



Article submitted to Royal Society 



T^X Paper 



2 



Madhusudhan VenkadesarJ and L. Mahadevan} 



of initial noise. Second, this sets the stage for understanding strategy and learning 
using a trial-to-trial iterative process of error estimation and correction. 



Not surprisingly, throwing accuracy has been studied by many others in var- 



ious 


contexts such as clinical applications (Timmann et al. 


1999 


Powls et al. 


1995 


, human evolution 


Woodetal. 2007 


Watson 


2001 


Westergaard et al. 


2000 


Chowdhary and Challis 


19991, sports (Gablonsky and Lang 


2005: f 


Smeets et al. 


2002 


Dupuyetal. 2000 Brancazio 1981 


Tan and Miller 


1981 OkuboandHub- 


bard 


2006), and human motor control (Hore et al. 


1996 Kudo et al. 


2000 


Cohen 


and Sternad, 2009). However, none of them address the choice of throwing style, or 



present an analysis of error propagation through projectile dynamics in the presence 
of an arm, which we consider here. 



In Section [2] we outline the various possible physical and biological influences 
on throwing strategy. Section [3] develops our model of throwing along with the so- 
lution for a perfect strike, while in Section [4] we use a linearization of the mapping 
between release parameters and projectile landing location to quantify error ampli- 
fication. In Section[5]we calculate the dependence of error amplification on throwing 
speed. Section [6] presents the a fully nonlinear calculations for propagating noise in 
release parameters with assumed distributions, and compares the linear and non- 
linear predictions using a specific numerical example. Section [7] analyzes the limit 
of vanishing arm length, i.e. the problem of shooting a projectile and calculatesthe 
optimal shooting angle as a function of the relative noise level in shooting angle 
versus speed. We conclude in Section [8] with a comparison of our findings with 
previous results, and some implications of our results. 

2. Factors that affect throwing style 

The choice of strategy for any motor task is influenced by a number of physical 
and biological factors. For throwing, the physics of projectile flight is affected by 
the projectile's shape, air drag, spin, wind, and so on, which in turn affects error 
propagation by flight dynamics. The decision of how to release a projectile is de- 
termined by a number of biological factors including muscle strength, skeletal joint 
limits, sensory feedback, memory, and learning. 

Here, we only study the case of windless, drag-free point-like projectile with no 
spin. The target can itself have many variations such as its shape, orientation, size, 
possibility of rebounds, a backboard (like in basketball), and so on. We consider 
two target geometries: (i) an upward facing, horizontal line-like target of prescribed 
small extent, and (ii) a circular target of prescribed small extent. 

The complex articulation of the human arm allows it to control four independent 
launch parameters in the plane of projectile flight, for example, the location and 
velocity of the projectile. Here we model the arm much more simply as a hinged 
rigid bar connecting the shoulder to the hand, so that it can only control two 
independent launch parameters. Although the neural controller can control any 
two independent launch parameters, say the speed and angle of release, or the 
speed and timing of release of the projectile, and so on, we use the angle the arm 
makes with the horizontal (4>), and the angular velocity (fi). We also disregard any 
strength differences between throwing styles. We note that prescribing an alternate 
set of variables is not completely equivalent to our choice, since the Jacobian of the 
transformation will necessarily play an important role in determining how errors 
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Figure 1. Model of the arm used to investigate how errors in throwing are propagated by 
projectile mechanics. The arm is modelled as a hinged rod of fixed length whose angle 
and angular velocity can be fully controlled. Using this 'super-human' arm model that 
has no limits on strength or coordination ability allows us to focus only on the projectile 
mechanics without consideration of ability to control the arm. R is the arm's length, cp and 
il are the angle measured from the +X axis and angular velocity of the arm, respectively. 
L and H specify the target coordinates, and g is gravitational acceleration. 

are amplified. This implies an important role for the choice of representation of the 
control variables for the motor task, although here we will limit ourselves to just a 
single convenient one given above. 

The variability and correlations in the release parameters (<p, f2) depend on 
the detailed properties of the neuromuscular system. For simplicity, here we use 
uncorrelated noise in <j) and tt. For the linear analysis we assume infinitesimally 
small noise in (f> and il and therefore need no further assumptions about the specific 
distribution. 




3. Model of throwing 

With the simplifications outlined above, we model the arm as a hinged bar of 
fixed length with independent control over the arm's release angle ((f)) and angular 
velocity (fi) as shown in Figure [T] For a drag-free point-like projectile moving in 
a uniform gravitational field, its trajectory (X(T),Y(T)) as a function of time T 
measured following release is obtained by solving the initial value problem, 

X(T) = 0, X(0) = -Rflsinct), X(0) = i?cos0 (3. fa) 

Y(T) = -g, Y(0) = Rttcosct>, Y(0) = Rsmc/) (3.1b) 

where, X and Y are measured relative to the arm's pivot, R is the arm length, 
and g is the acceleration due to gravity. We use the arm's length R and its natural 
frequency \J g/R to define the following dimensionless variables, 




Then the dimensionless equations for the projectile to land on a horizontally ori- 
ented, planar target (like a bin) located at a scaled height h = H/R, and distance 
I = L/R are, 

x(t, 4>, w) = cos <j) — tu sin <j) = I (3.3a) 
y(t, 4>, oj) = sin <p + tea cos <j) 1 2 = h (3.3b) 
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for t, one equation with three unknowns, and using 
- 1 < for the projectile to strike the upper face of 
the target, we obtain a solution surface th(<p,oj). By substituting t = ti l ((f>,uj) in 



By solving equation 3.3b 
the condition y(t) — lu cos <j> 



equation 3.3a| we also obtain an equation for the horizontal landing location of 
the projectile when it strikes the plane of the target. 



th(<j>,(jS) — u> cos (j) + \J to 2 cos 2 (f> — 2(h — sin <f>) 
Xh((j>,u) = cos <j> — w sin (f> cos (f> + yj uj 2 cos 2 (f> — 2(h 



— sin t 



(3.4a) 
(3.4b) 



Note the necessary condition that lj 2 > 2(h — s'm(f>)/ cos 2 <f>, i.e. a minimum launch 
velocity is needed to reach the target's plane when h > sin (f>. 

Then we find the conditions for exactly striking the target by solving equa- 
tions 



3.3 simultaneously, or by setting Xh( 



one parameter family of solutions 



I in equation 3.4b to obtain a 



cos < 



sin « 



2 

sin <f> 



(1 — lcos<j> — hsint 



(3.5) 



as shown, for example, in Figure [2t 



4. Quantification of error amplification 

When the release parameters (</>, u>) deviate from the one-dimensional solution curve 
((f), u)q (4>)), the projectile misses the target leading to an error Sxh = Xh(4>, ui) — l. To 
quantify the amplification of small launch errors in the throw, we linearize x^(</>, to) 
in the neighbourhood of the curve ((^,wo(^)) to obtain the amplification of 'input 
errors' (6<j> 6ui) T to the 'output /target error' Sxh- 



5u> 



where J P 



(A\ _ ( dx h dx h 
c\9) — I d4> du> 



(4.1a) 
(4.1b) 



The amplification of errors as a function of <f> is quantified by the largest singular 
value (A(</>)) of J orr (0) (The explicit expressions for J e „((j)) and A(</>) are given in 
Appendix A ). Because there is a one-dimensional curve (</>, wo(</>)) where Sxh = 0, 
J orr will be rank deficient, i.e. it has a zero singular value and an associated non- 
trivial null-space, namely, the tangent to the curve (4>,Uq ((/>)). 

By comparing the lowest order Taylor series expansion of (Sxh) 2 given by 



{5x h f 



5u> 



r I 

u err u err 



8ui 



Slu 



(4.2) 



with equation 4.1a , we see that the quantity 2 A 2 is equal to the largest eigenvalue 
of the Hessian (H eIT ) of the output variance (Sxi) and thus the maximal princi- 
pal curvature of the surface Sx 2 ((f), wq (0)) = shown in Figure [2]d, leading to a 
geometric interpretation of error amplification. 
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Figure 2. One-dimensional curve of parameters for an accurate strike, and error amplifi- 
cation, a. For a given target location (I = 1.5, h — —1.5) there is a curve of parameters 
((j>,u}o(4>)) for exactly striking the target. The portion with positive cj is the underarm 
throw according to our sign conventions and uj is negative for the overarm throw. There 
are some angles of release (</>) for which is it impossible to strike the target as seen by 
the non-overlapping region between the underarm and overarm throws, b. Deviations 
away from the curve leads to an error 5xh, and we plot the squared error. Local error 
amplification near ((f>,uo(4>)) is quantified by the maximal curvature that varies with <f>. 

To quantify the relative grading of strategies, we start with a simple numerical 
example shown in Figure [3ji which compares two throwing strategies — overarm 
(red) and underarm (green). Equal errors in release parameters amplify differently 
for different strategies, because A for the overarm throw is smaller than for the 
underarm throw for that target location. The envelope of trajectories, and hence 
landing locations of the projectile, was found by integrating forward a large number 
(100) trajectories starting from initial conditions in a box (0 ± 8(f), u> ± Suj). The 
nominal values for <f> and u> for the two throwing styles were chosen as the 'optimal' 
strategy in the sense we define below. 

Since the amplification of input error is directly proportional to A, a natural 
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metric to grade strategies for a given target (I, h) is the reciprocal of A, which we 
will henceforth refer to as 'throw accuracy', or simply 'accuracy', 

K^) = T7^» 0e{^:<5x h (0',woW)) = O} (4.3) 

In Figure [3J3 we show a polar plot with polar-angle tfi and polar-radius p(<f>), for 
two different target locations — one below and another above the arm's pivot. Recall 
that u>o(4>) > for underarm and wo(</>) < for overarm. For infinitesimal errors 
in release parameters, the best underarm throw is clearly superior for the target 
above the pivot, and the overarm strategy is better for the target below the pivot. 
If one were not a perfect planner, i.e. one does not choose the throwing angle cf> 
that corresponds to the maximum, the answer is no longer as clear. In Figure 
for example, the accuracy of the underarm throw in the right polar plot has a very 
narrow peak, i.e. small deviations in planning the angle of release from the optimum 
cause a drastic deterioration in accuracy. Therefore, although the underarm throw 
is superior for small planning errors, it is not robust to large planning errors. By 
accounting for both the peak and width of p(4>) we predict the fraction of overarm 
throws for a given target, and compared it against data in a separate publication 



(Vcnkadesan and Mahadevan In Review). 



All of the above calculations are for a flat, horizontal target, like a bin. How- 
ever, if the target had a different shape, then target error will no longer be the 
horizontal deviation Sx^- For example, consider a ball-like circular target, where 
the target error will be the distance of nearest approach (Figure B. Unlike a flat 
target (Figure j3pV) , the locus of nearest approach points for a circular target de- 
pends on the initial conditions (Figure |4| . These differences in error propagation 
naturally bear out as sensitivity of the optimal throwing strategy to the shape of 



the target (Venkadesan and Mahadevan In Review). 



5. Speed vs. accuracy 



We have so far presented accuracy p as a function of release angle <j> in the neigh- 
bourhood of the curve of accurate strike solutions u>o(<f>). In order to understand if 
there is a relation between speed and accuracy, we re-parameterize p as a function 
of uj by inverting the function ujq (</>) as shown in Figure [5][f] We see that slower 
throws are typically more accurate than faster ones in our model — the classic trade- 



off between speed and accuracy that is observed in human motor behaviour ( Fitts 



19541 |Meyer et al.| [19901 |Kerr and Langolf] [19771 |Harris and Wolpert| [19981 



Todorov 



2004 



2006 



Etnyre 1998). Such speed-accuracy trade-offs in biological systems 



has so far been explained as the result of signal-dependent and structured covari 
ance of input noise at the level of muscles ( Harris and Wolpert[ 1998 2006 Todorov 
2004). Our model however, has no strength limitations, and the input errors are 



additive, uncorrelated and independent of posture/velocity. Yet, we see the emer- 
gence of a speed-accuracy trade-off dictated only by projectile dynamics, and is 
thus a surprise. 

f We note that the function uj = ujq(</)) is not bijective (Figure|2^,) , and generally, the pre- image 
of every uj consists of four distinct values of c/> corresponding to a high/shallow, overarm/underarm 
throw (four curves at each location in Figure [5j). 
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Figure 3. Errors in releasing the projectile are amplified differently for different throwing 
strategies, depending on the target location and planning errors, a. This figure shows a 
comparison between a specific overarm (solid) and underarm (dashed) throw by the simple 
hinged arm. Throws where u)o(<f>) > are called 'underarm' and those with u>o{(f) < are 
called 'overarm'. For the same initial errors in angle 8(f) and angular velocity Soj of the 
arm at release, the output error 5xh for the overarm strategy is lesser than the underarm 
strategy. We quantify this error amplification using the non-zero singular value X(4>) of 
the map J erT : (54>,5lj) t h» 8x h . b. Polar plots of accuracy p(4>) = l/\(<j>) as a function 
of arm angle at release <j) f° r two different targets. If the choice of release angle is close to 
optimal, the underarm throw is clearly superior for the target above the pivot, and the 
overarm strategy is better for the target below the pivot. But, in right polar plot, a small 
(but non-infinitesimal) deviation from the optimum will cause a drastic deterioration in 
accuracy for the underarm throw, although it is superior in terms of its maximum. 



For every target location, there exists a value u> m i n (l,h), which is the smallest 
velocity needed to reach the target. For faster throwing velocities to > w m i n , there 
are four distinct throws that hit the target — overarm/underarm, shallow/ high 
throw. We find that the shallow overarm throw is most accurate for higher speeds, 
particularly for targets at or below the arm's pivot (equivalent to the shoulder). 
However, for slower speeds, the underarm throw is marginally more accurate than 
overarm for targets above the arm's pivot. 
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37°±6° 
-0.76 ±0.01 



Circular 
target 



: 280° ±6° 
u) = 1.10±0.01 



1.5 
-1.5 



Figure 4. Errors for a circular target depend on the angle of approach of the projectile 
as seen in the numerical example shown above. This is because, error is defined as the 
distance of nearest approach. 

6. Propagating distributions with non-infinitesimal variance 

We next relax the assumption of infinitesimally small input errors that was needed 
for the linearized analyses presented thus far. We introduce finite noise at the input 
side using a joint probability density function (PDF) /u <*,] associated with the noise 
in (j> and U). The horizontal distance x% at which the projectile lands is given by the 
non- invert ible, surjective function Xh(<f>,w), which transforms into the PDF 

f Xh associated with x%. For a fixed <f>, fu tU1 ]((j),u)o((j),l)) is transformed into f Xh {l) 
by the Jacobian of loq with respect to I. However, for a given target at distance 
I and height h, there is a curve of release parameters (<p,u)o(<p,l)) for accurately 
striking the target. Therefore, integrating over this solution curve and using 
for uj , we find f x as, 



3.5 



du>o(tj>' ,1) 



du>g((j)',l) 
dl 



d(j>' 
d(j)' dl 



— oo <3> 

Oluq (4>,l) esc 2 4>{4h sin (f> + 21 cos <f> + cos 20 — 3) 



01 



4V2(cot 4> + esc (f> — h — 3/2 



(6.1a) 
(6.1b) 



where $ = {</> : x^c/), u>o{<j), I)) = 1} is typically disjoint and depends on the location 
of the target as enumerated in |Appcndix B To ensure that f Xh is a probability 
density with total area equal to one, we normalize 6.1a| [f] 

To follow the implications of the fully nonlinear calculation for error amplifi- 
cation, we choose the simple example of throwing a projectile into a bin. For this 
example, we assume uncorrelated noise in <fi and lo, fu,u] — f<pfu, where f$ is a 
von Mises distribution. Because </> is a periodic variable, we use a von Mises distri- 
bution (von Mises 19181 iMardia and Jupp 20001, which is the circular analogue 



f There is 'leakage' because uj) is not a real number for some combinations (</>,lj) when 

the projectile fails to reach the plane of the target, and thus normalization becomes necessary. 
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Figure 5. Speed-accuracy trade-off for various locations of the target. Each sub-plot shows 
p(cjn) as a function of uiq for a different target location. One could hit a given target with 
many possible velocities by changing the angle of release. However, each launch velocity 
(and hence, each launch angle) has an associated error amplification quantified by p = 1/A. 
As seen in this figure, slower throws are typically more accurate than faster ones according 
to our model, much like the classic speed-accuracy trade-off that is seen in human motor 
behaviour. It is noteworthy that for some target locations, our model predicts that the 
most accurate throw is at speeds slight higher than ctj m i n the smallest possible speed to 
reach the target. This feature is more pronounced for nearer targets that are below the 
arm's pivot. Note that for uio > cj m i n , there are 4 different ways to hit the target, a shallow 
or a high throw and an underarm or overarm throw. This figure shows p(u>o) for each of 
these 4 throwing strategies by restricting our attention to uo £ [w m i n , 3.54]. For a human 
arm of 1 m length, cj = 3.54 corresponds to a launch speed of ~ ll.lm/s « 40km/h. We 
also find that when higher speed throws are necessary, a shallow overarm throw is more 
accurate than any other strategy. 

of a Gaussian distribution. The PDF of a von Mises distribution with mean [i and 
concentration k (1/k is analogous to the variance of a normal distribution) is given 



where Iq(k) is the modified Bessel function of order 0. For f u , we use a one-sided 
truncated Gaussian in order to restrict u to a single throwing style at a time, i.e. 
lj > or ui < for underarm and overarm, respectively. The PDF of a truncated 
Gaussian variable x with mean /x and standard deviation a that is truncated at a 
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Table 1. Parameters values used in nonlinear calculations for bin at I = 2, h 



Overarm 


Underarm 


/* 








joptim 
^rovcr 








.opt™ 

' under 




/ i optim\ 




1.08 


0.09 


-1.19 


0.18 


5.72 


0.09 


1.01 


0.18 


62° 


5° 


3.72 m/s 
13.4 km/hr 


0.56 m/s 
2 km/hr 


-33° 


5° 


3.17 m/s 
11.4 km/hr 


0.56 m/s 
2 km/hr 



and b is given by, 



g{x,fj,,a,a,b) 







if a < x < 6 
Otherwise 



(6.3) 



where erf is the Gauss error function. Therefore, the distributions we use for nu- 
merical examples are, 



t>') — v((f>', 0°_^ tlm , l/o - ^), (■) = over or under 



/ (a/) = I w o(^o?CT m ), cr w , -oo, 0) for overarm 
1 fl , (w / ,w (</>°^™),CT w ,0,oo) for underarm 



(6.4a) 
(6.4b) 

(6.4c) 



where a is the standard deviation, and 0°P tlm denotes the optimal throwing angle 
such that the linear amplification A(</> optlm ) is the minimum for the chosen throwing 
style ( |A 2| and Sections El). 

We now perform linear and nonlinear analyses for a specific numerical example, 
a bin that is 2 arm lengths away, and 1 arm length below the shoulder, i.e. / = 2, 
h = — 1. For this bin, we use the linear error amplification A to find the angle <f> and 
speed L) for the optimal overarm and underarm throws. The input distributions for 
the nonlinear calculation are centred at the optimal release parameters found from 
this linear analysis. The variance of the distributions are arbitrarily chosen, but 
are equal for both throwing styles. Table [T] lists the numerical values used in this 
example. Dimensional values are shown for a one meter long arm. 

Our numerical solutions show that the underarm throw (green solid curve in 
Figure [6]) amplifies errors more than the overarm throw (red solid curve in Fig- 
ure [6]), as predicted by the linear amplification A. Quantitatively as well, the ratio 
of overarm to underarm linear amplification is similar to the ratio of a, the stan- 
dard deviation of the output distribution f Xh (see table [2]). However, the nonlinear 
calculation calculation shows that the output distribution for an underarm throw is 
skewed (Figure [6j table [2J), something that cannot be found from a linear analysis. 
The implication of this skewed output is that the underarm throw could be per- 
ceived as a weaker throw because the mode of the distribution where the projectile 
is most likely to land, undershoots the target by 4% in this example. In contrast, 
the overarm throw undershoots the target by only 1%. We speculate that this un- 
dershoot, solely due to projectile dynamics, may partially underlie the common 
perception that an underarm throw is 'weaker' than an overarm throw. 
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Table 2. Comparison of linear and nonlinear calculations for bin at I = 2, h = — 1. 



Linear 


Nonlinear 


Aover 


0" over 


Mean landing location 


Most likely landing location 


^midcr 


u under 


Overarm 


Underarm 


Overarm 


Underarm 


0.69 


0.56 


1.99 


2.01 


1.98 
1% undershoot 


1.92 
4% undershoot 



(Linear: A(0°^ m ) = 1.35, \(<PZdD = 1.96; Nonlinear: cr ovcr = 0.19, cr undcr = 0.34.) 




1 2 3 



/ 

Figure 6. Propagating initial distributions with equal variance for hitting a target at I — 2, 
h = — 1. The underarm throw amplifies errors more than the overarm throw, in agreement 
with the linear analysis. For comparison, we superimpose a normal distribution (dashed 
curves) with mean and variance equal to the propagated distribution (solid curves), to- 
gether with the mode (solid vertical line) and mean (dashed vertical line). The overarm 
throw leads to a nearly normal distribution where the projectile lands with mode (1.98) and 
mean (1.99) nearly equal to each other. The underarm throw in contrast leads to a skewed 
distribution where the most likely outcome is an undershoot as seen from the mode (1.92) 
although the mean (2.01) is nearly at the target (1 = 2). The input distributions are centred 
at (0°£r m = 1-082, o;o(0C?er m ) = -1-186) for the overarm throw and at (<#£*j™ = 5 - 716 > 
w ° (^und™) = 1-008) for the underarm throw, and we set standard deviations as = 0.087, 
o~u — 0.177 for both throws. For a one metre long arm, these values expressed as mean ± 
standard deviation are, overarm: <f> = 62 ± 5°, \v\ — 3.72 ± 0.56m/s (13.4 ± 2km/hr) and 
underarm: </> = -33 ±5° ,\v\ = 3.17 ± 0.56 m/s (11.4 ± 2 km/hr). 

7. Shooting: Zero arm length 

Our considerations so far use the length of the arm R to set the length scale in 
the problem. We now look at the limit when this length scale becomes vanishingly 
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small or for throwing to faraway targets where L, H ^> R. In this, the artillery 
limit, the problem becomes one of shooting a projectile at an angle 9 with a linear 
velocity V. Naturally therefore, we only consider targets at the same height as the 
origin from where the projectile is launched. The ratio of variability in V to the 
variability in 9, introduces a velocity scale k = Sv/69, where the quantities 5v and 
89 are appropriate measures of variability. By expressing distances in units of k 2 / g 
and time in units of k/g, we obtain the scaled trajectory equations for distance x, 
height y, as a function of time t, launch angle 9 and launch velocity v, 



x(t) = tv cos 9 

y(t) = tvsm9 - \t 2 



(7.1a) 
(7.1b) 



For exactly striking the target at a distance I with the launch parameters 6,v, the 
projectile lands at xq(9,v) = I = gL/k. This gives the one-dimensional curve of 
launch parameters (9q,vq(9q)): 



vo(0) 



I 



sin26» 



(7.2) 



Performing a linearized error analysis in the neighbourhood of this curve (#0, vq(9q)) 
yields the Jacobian J that maps small errors 69 and 5v to small errors Sxq- 



e=e ,v=v (e ) 



= [2vl cos 26* 2v sin 29) = 2 i 



sin 26 



1 VI sin 26>) 
(7.3a) 



The shooting angle most robust to small input noise corresponds to the 9 where 
the only positive singular value is minimized. Minimizing the square of the singular 
value 

C(0) = J T J= l ^°i 2e +lsm29 



yields 



d0 

> cos 29* 



sin 2 29 



O^Al cos 26»* - 81 



, cos 26»* 
sin 3 26 1 * 



= 



0, or, for 0< I < -, sin 20* 
2 



(7.4) 

(7.5a) 
(7.5b) 



We therefore find that for I > 1 /2, there is only one extremum at 9* — ir/4, but three 
extrema co-exist for < I < 1/2 at 9* = 7r/4, arcsin \/2l/2, and tt/2 — arcsin v2Z/2. 
To determine if the solution |7.5b] is a minimum or maximum, we evaluate 



ft 
d9 2 



-16Zsin 29 



321 2 + 6il 2 cos 2 29 



ft 
d9 2 



sin 4 26* 

j 161(21-1) when i 

1 24(22)3(1 - (2/) I) when sin 29* = \/2l, < I < \ 



f > < I 



(7.6a) 
(7.6b) 



From the solutions in |7.6b| we see that 9* = tt/4 is a local minimum for I > 1/2, a 
local maximum for 2 < 1/2, and there is a pitchfork bifurcation at I = 1/2. We also 
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Figure 7. Bifurcation diagram for best shooting strategy based on a linearization in the 
neighborhood of all strategies that lead to an accurate strike. Recall that I = gL/k 2 is the 
dimensionless distance to the target and ^optimal is the angle to shoot the target so that 
launching errors are amplified the least amount. 



find that for I < 1/2, when 0* = 7r/4 becomes a local maximum, there are two new 
branches that are a local minimum, with a cubic dependence on I near I = 1/2. 

The result that the optimal shooting angle (0*) for targets with I > 1/2 is 7r/4 
is well known, including cases where there is air drag (Gablonsky and Lang 2005 1 . 
However these past results however do not scale the equations like we do using 
the relative amount of noise in shooting velocity compared to shooting angle, and 
therefore do not identify the pitchfork bifurcation. Thus = ir/4 becomes the worst 
possible choice for shooting when I < 1/2, and there exist two symmetric branches 
of optimal shooting angles. 

We attempt to provide physical intuition to this surprising result for the problem 
of shooting a drag-free projectile. In the limit of the relative noise in shooting 
velocity being much greater than that in the shooting angle, i.e. k 3> 1 =>■ I ~ 0, 
the best shooting angle corresponds to zero sensitivity to fluctuations in velocity, 
namely a high velocity shot with 9 = or = n/2. At the other extreme where 
fc w =^ i > 1, i.e. when the noise in shooting angle is much larger than in 
shooting velocity, 7r/4 is the optimal angle because that corresponds to a local 
maximum for projectile range, and hence has zero sensitivity to fluctuations in 
shooting angle. There is a bifurcation, and not a smooth continuous change between 
these two extreme cases because the problem is symmetric about 7r/4 and the 
optimal shooting angle depends continuously on I. The specific bifurcation point 
is determined by a trade off between the two sensitivities of the shooting range to 
shooting angle and velocity, respectively. For small perturbations e near 9 = it /A, 
the sensitivity to shooting angle is 16l 2 e 2 + 0(e 3 ) and the sensitivity to shooting 
velocity is — 8Ze 2 + 41 + 0(e 3 ). Therefore, in the neighbourhood of = tt/4 and for 
I < 1/2, the sensitivity to shooting velocity overpowers the sensitivity to shooting 
angle. 
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8. Discussion 



Our calculations of optimal throwing strategies are based on a consideration of how 
infinitesimal noise in release parameters is propagated by the parabolic trajectory 
of a projectile. Past studies have found that greater error is permitted in throwing 
angle, than speed ( Gablonsky and Lang| 2005 Brancazio 1981). However, there 
is an inherent problem in how to compare two quantities that are dimensionally 
different. Past researchers either treat errors in release angle and velocity separately 



(Dupuy et al. 20001, or compare percentage errors (Gablonsky and Lang 2005 



Brancazio 1981 ), quantities whose scaling is affected by the magnitude of throwing 



angle and velocity, how angle is measured and the relative weight of these measures. 
By introducing a natural length scale, the arm length, and making our formulation 
dimensionless, we can directly compare sensitivity to errors in angle and speed, thus 
avoiding multiobjective optimization that uses a heuristic weighting parameter. In 
our minimal model, we ignored sensory feedback from the arm, covariance structure 
of the noise in input parameters, dependence of noise on projectile mass, throwing 
speed, and so on. Nevertheless, we find a dependence of optimal strategy on target 
geometry, and a trade-off between speed and accuracy, which has until now been 
ascribed to signal dependent noise in muscles (Fitts 1954). Our results agree with 
the prevalent observation in the literature that the slowest throw is most accurate 



(Gablonsky and Lang 2005 Brancazio 1981 Tan and Miller 1981 Okubo and 



Hubbard 2006). A minor modification that accounts for the distribution of angles 



associated with planning, allows our theory to make some testable predictions for 



the fraction of overarm throws used by humans in an experiment ( Venkadesan and 



Mahadevan In Review). 



We have also identified a non-trivial null-space of the mapping between errors 
in release parameters and errors in striking the target. This naturally suggests 
an additional control strategy to maximize accuracy of throwing: by tuning the 
covariance of input errors in angle and angular velocity to co-vary along the null- 
space of J crr for small noise, or along the curve (</>o, wo(</>o)) for large noise, one 
could mitigate output errors. In fact, there is some evidence that humans control 
the covariance of the projectile's launch angle and linear velocity in such a manner 
(Dupuy et al. 2000 Kudo ct al. 2000 ). More generally, there is substantial evidence 
at the behavioural and neurophysiological levels that humans are adept at tuning 



covariance of noise across multiple joints / muscles (Todorov 2004 Valero-Cuevas 



et al. 2009). However, because it is likely impossible to perfectly covary input 



errors, our calculations of maximal error amplification should predict strategies for 
accurate throwing even in the presence of structured noise. To see this, we note 
that geometrically, the effect of structured noise amounts to asking how an ellipse 
gets amplified by J orr . Even if the long axis of the ellipse is oriented parallel to the 
null-space of J er r, so long as it has any non-zero length along its minor axis, the 
non-zero singular value of J err specifies how input errors get amplified. 

Finally an additional conclusion is that the most accurate throw is close to, but 
slightly faster than the slowest throw, an observation that has also been experi- 



mentally reported (Freeston et al. 2007). Moreover, we find that at higher speeds, 



there is a clear advantage for the overarm style. This is particularly salient in light 
of hunting and human evolution. Humans are significantly better (faster and more 
accurate) at overarm throwing than any other member of the animal kingdom 
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Wood et al. 


2007 


Watson 


2001 



Westergaard et al. 20001. Given that hunting 



possibly played a role in the evolution of human form, could the physics of error 
propagation in throwing provide a partial explanation for the selective pressure for 
superior overarm throwing skills in humans ? 
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Appendix A. Linearization for throwing with an arm 

When throwing with an arm (Sections [3j|5]), the Jacobian of Xh{4>,<^>) in the neigh- 
bourhood of the curve (</>,wq(0)), and its non-zero singular value is given by, 



M) = ( 



<fi,uo(<t>) 



dx h 
d(/> 

dxh 



Ahl cos <j> - A(h + I cot <j>) + esc 4> ( (2l 2 - l) cos 2(j> + 3) 
^ ^ Ah sin 4> + 21 cos </> + cos 20 - 3 

Ay/2 sin 2 <p(csc tj> — I cot <j) — h)i 



(Ala) 
(Alb) 
(Ale) 



The non-zero singular value A(0) of J e rr(0) is given by, 



A 2 W 



dx h 



4>^o (4>)/ 



dxh 
da; 



(A 2) 



<t>,uo (<t>)/ 



Appendix B. Feasible arm angles 
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Figure 8. The individual cases for identifying the feasible arm angles in order to be able 
to strike a planar target, a. Two qualitatively different shallow shots to hit the target 
depending on whether the target is above or below the point of release, b. Enumeration 
of different regions of target location where domains of <f> are calculated for Xh{§-, c^o(0)) 
to be a real number. 



To calculate the arm angles where it is possible to hit the target, we first observe 
that there are two qualitatively different shallow shots depending on the height of 
the target relative to the release point (Figure |8^i) . We refer to the straight-line 
shot with infinite velocity by (jf^ ^ and the curved shallow shot which is a parabola 
passing through the target with its apex at the target, and tangential to the unit 
circle by <jf, -.. The subscript will denote the throwing style ('u' or 'o'). These angles 
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are give by, 



& ,,. 3 est . „ , 9 a ,,. a 



c(a) = (1 + 0tan 4 - - 4/itan d - +6tan^ - - Ah tan - + (1 - Z) (B la) 

max a where a is a real root of c(a)if I < 1 , , 

€ = < . ., , ^ 1 (B lb) 
I mm a it ( > 1 

min a if i < 1 ,„ 

0u= (Blc) 
I max a it « > 1 

<ft s — arctan 2(h, I) + arccos - (B Id) 

</>u = arctan 2(h, I) — arccos - (B le) 

arctan 2(h, I) = arg(7 + ih) (B If) 

d = y/P + h 2 (B lg) 



B la has only two real roots in the interval [0, 2ir]. The domains of feasible 4> so 



that Xh(4>,u}o(<p)) is a re &l number are given by, 



1 


Z < > 1 


$ = 


(O,0S]U(7r,#] 




2 


I < -i,\h\ < l 


$ = 


(O,0 S U )U(7T,^] 




3 


I < -l,h < -1 


$ = 


M)U(7r,#) 




4 


l>l,h>l 


$ = 


[^,27r)U[^,7r) 




5 


l>l,\h\<l 


$ = 


[^,27r)U(^ ,7r) 




6 


l>l,h<-l 


$ = 


«,27r)U(^,7r) 




7 


\l\<l,h>l 


$ = 






8 


\l\ <l,h<-l 


$ = 


((O>]U(0 S U ,2^- 


■0]) U ([-0,tt) U [2tt 


9 


-1 < I <0,0 < h <1 


$ = 


(O,0 S U )U[0» 




10 


-1<1 <0,-l<h<0 


$ = 


((0,V]U«,2tt- 


■0]) U ([V>,7r) U [2tt 


11 


0<l<l,0<h<l 


$ = 


(o,^]u(Ct) 




12 


0<l<l,-l<h<0 


$ = 


((0,V]UK,2tt- 


0]) U ([V»,7r) U [2tt 


13 


I 2 + h 2 < 1 


$ = 


(0,^]U[^7r) 





where t/> = arccos(Z), and the domains are listed in the format (underarm) U (over- 
arm) . 
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